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Abstract 

We show that if we require positive definite probabihties, then frequently cited 
results of Prato, of Lenzi and Hilhorst on the nonextensive equilibrium statistical 
mechanics of gases are valid only for a limited range of the Tsallis parameter, q. 
We determine the range of validity of the Hilhorst and the Lenzi formulae. We then 
use various integral representations of the Gamma function to derive new formulae 
for one-dimensional gases whose range of validity is wider than that of the Hilhorst 
and the Prato formulae. We then apply the new formulae to the classical ideal gas 
and the Tonks gas. 
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1 Introduction 



In 1988, Tsallis postulated a new entropy function [1]. He was, in part, mo- 
tivated by the difficulties of the standard thermodynamics with nonextensive 
systems. Since that time, his methods have found various applications. (See 
references [2] to [8].) 

One of the main difficulties in using the Tsallis formalism is the proper incor- 
poration of the constraints. The calculation of thermodynamical quantities in 
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this generalized scheme is more involved than in the usual statistics, frequently 
making necessary the recourse to numerical approaches [13]. 



Analytic approaches, however, have great value because because these may 
suggest new approximation schemes where numerical methods might fail. An 
example of these analytic methods are integral transform methods. 

It has been claimed that the Tsallis partition function is an integral trans- 
form of the corresponding q=l partition function. TsaUis[12], Prato[13], and 
Lenzi[14] have developed various integral transform formulae. 

Intergral transform formulae [13] [10] [12] have been used in the past to study 
the TsaUis statistical mechanics of the classical ideal gas [13], black body ra- 
diation [2], and two-dimensional self-gravitating systems[10], among others. 
If Hilhorst-type formulae are invalid for some values of q, then these systems 
need to be restudied for these values of q. 

Our task is to reexamine the derivation of the integral transform formulae, and 
to examine their respective domains of validity. In the following discussion, we 
concentrate on the canonical ensemble. In particular, we use the unnormalized 
constraints formalism [16]. 

1.1 The Hilhorst Formula 

The Tsallis partition function for the canonical ensemble is 



Note that this is a restricted sum; the sum is only over all states for which 
I3{q — l)Hi > 0. The i^j's are the eigenvalues of the Hamiltonian. 

It is now convenient to introduce the Gamma function so that the power law 
form may be written as an integral of an exponential. The representation of 
the Gamma function to be used is given by 



(1) 



oo 




(2) 







where Re(x) > 0, Re(y) > 0.. 



Prom (2), it can easily be seen that, for Re(x) > 0, Re(y) > 0. 
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Setting X = 1 + P{q — l)Hi, and y = ^^j, the expression for pi may be written 
as 



oo 

J ^fEf e-*[i+/3(9-i)^^,ldi (4) 



where Re(l + - l)Hi) > 0, q > 1. 

Previous workers then went on to substitute (4) into the partition function to 
obtain 

oo 

^. = Ef7Io/^^e-WM^^W (5) 

i \q-l) Q 
oo 

= -^Jt'^^^e-^l^^^(^-^^^Mt (6) 

yq-l) i 



oo 

XT / t^e-*i:e-*[^(^-^)^^ldi (7) 



r(^) 



Previous workers then went on to write 



oo 

I t^^e-'Z^{t{q - l)(5)&t (8) 



Za = - 

Vg-l^ 



This is called the Hilhorst formula. 



1.2 The Lenzi formula 



In this section, the derivation of the Lenzi formula[10],[14] is reviewed. 

To obtain the Lenzi formula, consider the following representation of the 
Gamma function [18]: 
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exp{ab)b 



l-z °? r.itb 



[ , ^ dt (9) 
J a + itV ^ ' 



Viz) 27r J ia + ity 

—oo 

with a > 0,6 > 0, Re{z)> 0, and -7r/2 < arg(a + it) < 7r/2. 

Setting a — 1, b — 1 + P{q — 1)H, and z — -^j^ + 1, it may be seen that 

oo 

[1 + - l)Hi]^^ = I K,(i)e(-(^+")(^-^)^^)di (10) 

—oo 

where 



_ r[(2 - g)/(i - g)] exp(i + it) 

27r (i + ii)(2-.)/(i-9) ^-^'^ 

From the form of the partition function (16), and from (10) the TsaUis parti- 
tion function was written as 



Z.^Y. j i^,(i)e(-(^+")(^-^)^)^di (12) 

* — oo 
oo 

= j X,(i)^e(-(i+")(i-«)^)^di (13) 

— oo * 

The expression ^^e*^"^''^'^^*^^"^"'^^''^^ was then identified in previous works as 



Zi((l + it){l - q)P) = Y^e^-(^+''^i^-<i)m (14) 

i 

Using (14), the q < 1 partition function was written as 

oo 

Zg= J Kg{t)e-'Z,{{l + it){l-q)P)dt (15) 

— oo 

This expression is called the Lenzi formula. 

The Prato formula [13] was obtained in a manner similar to the way the Lenzi 
formula was obtained. The only difference is in the expression for the Gamma 
function. The Gamma function Prato used was expressed as a contour integral, 
and it may be shown (by a suitable deformation of the contour [19]) that the 
Prato and Lenzi formulae are equivalent. 
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2 How one can go wrong with integral transform methods 

We now reexamine the range of validity of the various integral transform meth- 
ods. As we have stated, integral transform formulae are important because 
they allow us, when these formulae are valid, to express the Tsallis partition 
function of a given system in terms of the g = 1 or Boltzmann-Gibbs partition 
function. Unfortunately, as we now show, integral transform methods are valid 
only under some restrictive conditions. 

These restrictive conditions have their roots in the positive-definite probabil- 
ity requirement. In deriving the Hilhorst formula, the sum J2i c"*'^*^''^^)^' was 
identified as the q = 1 partition function, with (3 replaced by t(3{q — 1). Only 
when this identification is valid may we deduce (8) from (7) 

This identification will not be valid for all cases because the sum in (2) is a 
restricted sum. That is, we only include states for which [1 -|- P{q — 1)H] > 0. 
On the other hand, for the q = 1 partition function, we have an unrestricted 
sum. Thus, only when the unrestricted sum is equal to the restricted sum will 
the integral transform formulae be valid. 

The Hilhorst formula will be vahd when q > 1 and all the eigenvalues of the 
Hamiltonian satisfy the inequality Hi > —-p^^zj;^- The set of positive definite 
Hamiltonians will obviously satisfy this requirement. Even if the Hamiltonian 
is not positive definite, there will be a range of values of g > 1 and (3 for which 
the Hilhorst formula will be valid. The Hilhorst formula will be vahd when the 
greatest lower bound of the set of eigenvalues of the Hamiltonian is greater 
than or equal to 

In the derivation of the Lenzi formula, the sum ^^e*^"*^^^^*^*^^^'''"^''-^ was identi- 
fied as the q = 1 partition function where (3 is replaced by (1 + it)(l — 
But this sum is a restricted sum; we must sum only over all states for which 
1 + I3{q — 1)H >{). Only when the restricted sum is equal to the unrestricted 
sum will the formulae of Prato and of Lenzi be valid. 

The Lenzi and the Prato formula will obviously be valid when g < 1 and the 
Hamiltonian is negative definite. In fact, Lenzi and the Prato formulae will be 
valid if g < 1 and the least upper bound of the eigenvalues of the Hamiltonian 
is less than or equal to j(YZq^- 

In general, if the Hamiltonian has both positive and negative values, we may 
not use the Hilhorst, the Prato nor the Lenzi integral transform formulae. To 
be able to apply integral transform methods, we need to derive formulae that 
work under less restrictive conditions. 
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3 Integral Transform Formulae for One-Dimensional Gases 

The derivation of more general formulae for one-dimensional gases is done 
here. Along the way, we will show an interesting feature of TsaUis statistical 

mechanics: for some values of /3, inaccessible regions of phase space may begin 
to appear. These inaccessible regions appear because of the positive probability 
requirement. The appearance of these inacessible regions in phase space is the 
reason for the failure of integral transform methods for some values of q. 
Thus, previous workers who used the integral transform methods neglected, at 
some point in their derivation, to take the positive probability requirement into 
account. This means their results need to be reexamined. 

In the following discussion, we concentrate on the canonical ensemble. In par- 
ticular, we use the unnormalized constraints formalism [16]. 



3.1 The q <1 case 

The TsaUis partition function for the canonical ensemble is 



The prime in the summation means we have a restricted sum; we sum only 
over all states for which 1 -|- I3{q — l)Hi > 0. The if^'s are the eigenvalues of 
the Hamiltonian. 

It must be noted that in previous work, the restriction of positive definite 
probabilities was not built into the notation. This ambiguity in the notation 
led to conclusions which are difficult to justify physically because the positive 
probability requirement was forgotten. 

Now suppose that we have a system of N identical particles in 1-dimension, 
with a Hamiltonian of the form 



Let us assume that the particles are confined in a box of length L. Then the 
partition function becomes 




(16) 



H^^^ + U{x^,...,XN) 



(17) 



i=l 
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= ^ / da:i... I d.x^ J dpi...dp^[l + (3{q - 1)H]'/^'-'^^ (18) 



again, the integration is only over states for which 1 + f3{q — 1)H > 0. Now 
suppose that the potential and j3 have values within the box so that 1 — 
l)U{xi, ...,xn) > 0. Then 



1=1 



n 1/(1-9) 



/3{l-q)El^t 
l-(3{q-l)U 



1/(1-9) 



(19) 



We now use spherical coordinates in momentum space. We set Y^f — P^, 
and note that the volume of an N-shell is dV{P) — Im/2)\ ^ ^-P- 

We may introduce the momentum cut-off explicitly in the following manner: 
Let 



PAP) 



2m[l - /3(1 - g)t/] 



(20) 



From the positive probability requirement, we wish to sum over all states for 
which P < Pf{(3). Thus we may write 



L L 

^'^m^ I ^""^-l '^''^ / dPi-dpiv[l + /3(g-l)//]'/('-^) 



L L 



ll/(l-9) 



I 



X / dP 



{N/2)\ 



1-9 



X 



Jdx,...l dx^(l-/3(l-g)C/)VM 

/ 



{N/2)\ 







(21) 
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Note that the integral over the momenta can be evaluated using the Beta 
function [19] 



± 

B{p, q) = '2j x^P-\l-xY 



-1 _ r(p)r(g) 
r(p + g) 



where Re(p) > —1, and Re(g) > —1. 

We then get, using the change of integration variables P — Pf{/3)x 



(22) 



/ 



dPP 



N-l 



1 - 



1-9 



K 1 



dx{l — x^) 1-9 x 



N-l 



= 2 



(23) 
(24) 



Substituting the expression in (24) into the expression for the partition func- 
tion (21), we obtain 



Z, = ^ / dxi... / da;^(l - P{1 - q)U{x,, a;^))^/2+V(i-.) 



x2 



2m \^/^„/iV 



(7V/2)! 
(25) 



Using the expression for the Beta function in terms of Gamma functions, and 
using the properties of the Gamma function, we may simplify (25) 



2m 



N/2 



47r^/2r(l/(l - g) + 1) 



V/?(l -q)J h^Nir{{l/{l - q)) + {N/2) + 1) 

L L 

X j dx,...J dx4l-P{l-q)U{xi,...,xr,)f/'+'/^'-'^^ 



(26) 



Note that this expression is only valid when 1 — f3{l — q)U >0 over the region 
of integration. If we wish to consider ranges of (3, or potentials U for which 
this restriction does not hold, we get the formula 
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^P{1 - q)J h^mr{{l/{l - q)) + {N/2) + 1) 

L L 


x^(l-/?(l-g)[/(,Ti,...,x,v)) (27) 
where the step function 9{x) is defined as 

9{x) = 1 when a; > 

= when x <0 (28) 

We now introduce the Gamma function because we wish to write a power 
law form as an integral of an exponential. The relevant representation of the 
Gamma function was used in the papers of Lenzi [14], [10], and is given by [18]: 

1 exp{ab)b^-' 7 e^*^ 



r(^) 2n J (a + it) 

— oo 

Prom (29), we have 

.1 T{z)exp{ab) J e^**" 



/ — -dt. where 6 > (29) 
J a + itY ^ ^ 



h 



-dt where 6 > (30) 

(a + itY ^ ' 

—oo 

Setting a = 1, 6 = 1 — /3(1 — q)U, and z = ^ + iZ^ + 1) we have 

r(f + ^ + 1) 7 i+it -(i+imi-q)u 

—oo 

Substituting (31) into (27), we get 



^ _f 2m 47r^/2r(l/(l - g) + 1) 



'1 



\(3(1 - q)J h^NW({l/(l - q)) + (N/2) + 1) 

f } r(f + ^ + l)y? Ql+itQ-{l+imi-q)U 
X / dXi... / dXN / ; r dt 

J ^ J 27r J 1 + ity 

-oo ^ ' 

xe{l-(3{l-q)U{xi,...,XN)) (32) 
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Interchanging the order of integration, we get 



/ 2m \^/' 47r^/^r(l/(l-g) + l) 
X / dt — ^ 

L L 

X J dxi... J dx;ve-('+'*)^('-^)^^(l - - q)U) (33) 



The presence of the step function in (33) tells us that if /? is large enough, 
then there ought to be inaccessible particle configurations. 

Note that this expression may be further simplified if the Hamihonian or the 
parameter f3 is such that the condition 1 — /?(1 — g)?7>0is satisfied for every 
point in configuration space. If this condition is satisfied, then the expression 
for the q < 1 partition function reduces to 



/ 2m \ ^^Z' 47r^/^r(l/(l - g) + 1) 
\P{l-q) J Wn\2^ 

°° l+it ^ ^ 

X / dt ^^^^ / dxi.../ dx^e-(^+^*)/^(^-^)^ (34) 

-oo (1 + it) 2 1 « Q 

/ 2,n \ 4^^^^^r(l/(l -q) + l) 
[p{l-q)) h^N\27r 

oo 

X / %' 1 ^l ^con/((l + ^t)/?(l-g)) (35) 

_l (l + it)T + Tr5+l 



Let us now discuss this very interesting result. First, the presence of the step 
function in (33) tells us that if (3 is large enough, then there ought to be 
inaccessible particle configurations. That is, there exist points in phase space 
that become inaccessible if we lower the temperature. This is in contrast to 
the usual Gibbsian statistics, where lowering the temperature decreases, but 
does not totally wipe out, the possibility of the system being in high energy 
states. 

Second, when 1 ~ f3{l — q)U >0 for every point in configuration space, then 
there is an exact relationship between the Tsallis partition function and the 
Gibbs configuration partition function. So, for potential energy functions that 
are bounded from above, there is a temperature range for which (35) is valid. 
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3.2 The q> 1 case 



We now consider the q > 1 case. Again we pay careful attention to the restric- 
tion of positive definite probabihties. The derivation is similar but not totally 
identical to the derivation of the Hilhorst formula [12] 

Again, consider the partition function (18) 



1 

Z, = I dxi... I dx^ I dpi...dp^[l + P{q - 1)//] VM 



(36) 



We must integrate over the region of phase space for which 1+ f3[q — 1)H > 0. 
Again, we introduce spherical coordinates so that we may write 



1 ^ L 



1/(1-9) 









^N/2]ypN- 



xeiu- 



{N/2)\ 



/3(i-g) 

1- P{l-q)U2m 



(37) 



and 



1 ^ ^ 
= ^ / dx^... I dx^il - /3{1 - q)U{xu ...^x^))'/^'"^^ 



X 



/ 

/ 2m[l-g(l~ 

V — 



dP 



^N/2jypN- 

W2)! 



/3(i-g) 
l-/3(l-g)C/2m 



1-9 



./3(l-g) 



(38) 



The full partition function is given by 



(39) 



Let us now examine our results. The first term in the expression (39) for the 
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full partition function reduces to the Hilhorst formula when the Hamiltonian 
is positive definite. 



The second term represents the correction to the Hilhorst formula when the 
Hamiltonian is not positive definite. Note that the presence of the step function 
tells us that whenever (^gp^^ ~ ^) > fo'^ ^ given point in configuration 
space, then there is a lower limit to the allowed kinetic energies. 

There is a problem with these expressions that we were unable to address: the 
convergence of the integrals in expressions (37) and(38). If the integrals in (37) 
and in (38)do not converge for a given value of q, then there might be reason 
to suppose that the probability distributions obtained are unnormalizable. 



4 Simple Applications 

4.1 The Ideal Gas 

The ideal gas is the simplest classical gas. 

We now reconsider the ideal gas within the unnormalized constraints formal- 
ism. Since the Hilhorst formula is valid for the ideal gas in the q > 1 regime 
and the since the q > 1 partition function has already been worked out cor- 
rectly [12], we work out the one-dimensional, q < 1 case. Since the potential 
is always zero, we may use the results of the previous section to immediately 
write the partition function as 



^ ( 2m \ 47r^/^L^r(l/(l - g) + 1) 

' W-q)J h^N\T{{l/{l-q)) + {N/2) + l) ^ ' 



On the other hand, Prato's formula[13] leads to the result 



2m \ ""'^ 27r^/2L^r(l/(l - g) + 1 



Ul-?)J h^N\T{{l/{l-q)) + {N/2) + l) ^^^^ 

Evidently, Z^"^"-^" < Zg. This is understandable because negative probability 
states are included in the expression for the g < 1 partition function. 



12 



4-2 The Tonks gas 



The Tonks gas [17] is a one-dimensional gas with a hard-core potential energy 
function. If the particles are a distance greater than a, then they act like free 

particles. However, they cannot be nearer to each other than the distance a be- 
cause the potential energy becomes infinite. In equations, the Tonks potential 
between two particles, (called the ith and jth particles) is given by 



Vronksixi - Xj) =0 \xi-Xj\ > a 

— oo \ Xi — Xj \ < a (42) 



It has been shown that the Tonks gas configuration partition function, in the 
Boltzmann-Gibbs treatment, is of the form 



r^,, _ {L-{N- l)ar 

^conf - ]yj ■ 



Note that the Hamiltonian of the Tonks gas is positive definite. Thus, the 
Hilhorst formula may be used to evaluate the q > 1 partition function. Sub- 
stituting the form of the Tonks gas Boltzmann-gibbs partition function into 
the expression 



oo 

I if^e-%(t(?-l)/3)dt, (44) 



r(^) 




the following is readily obtained 



= mr{^,) [m^)) ''^1 



5 Conclusion 



The positive probability requirement (or the need for cut-off conditions) is 
a feature of Tsallis statistics that is not found in the usual statistics. It is 
this requirement that makes Tsallis statistics interesting and computationally 
nontrivial. 
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We chose to work within the unnormahsed constraints formahsm purely for 
mathematical convenience. There are, however, other Tsallis formalisms [16] 
aside from the unnormahsed constraints formalism. There is therefore a need 
for a reinvestigation of these other formahsms in the hght of the positive 
probability requirement. 

We have shown, within the unnormahsed constraints formalism, how integral 
transform methods may fail for some values of the parameter q. The failure 
of integral transform method emerges from the appearance of innaccessible 
phase space configurations; these innaccesible phase space configurations in 
turn arise from the positive probability requirement. 

We then worked out expressions for the q < 1 and q > 1 Tsalhs partition 
function for one-dimensional gases made of identical particles that agree with 
the positive probability requirement. We were able to apply the formulae we 
obtained to the case of the Ideal Gas and the Tonks gas. 
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